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Abstract. Problem of cosmological singularity of general relativity theory 
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j>^! is discussed. The possible resolution of this problem in the framework of 

a ■ 

inflationary cosmology is proposed. Physical conditions leading to bounc- 
ing inflationary solutions in the frame of general relativity theory and 

gauge theories of gravitation are compared. It is shown that gauge theo- 
CN ■ 

ries of gravitation allow to build regular inflationary cosmological models 
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of closed, open and flat type with dominating ultrarelativistic matter at 
a bounce. 
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1 Problem of cosmological singularity and vacuum 

cr 

gravitational repulsion effect 



Problem of cosmological singularity (PCS) is one of the most principal problems of general 
relativity theory (GR). According to Hawking- Penrose theorems the appearance of cosmo- 
logical singularity in cosmological solutions of GR is inevitable, if gravitating matter satisfies 
so-called energy dominance conditions [1, 2]. There were many attempts to resolve PCS. If 
the Planckian (pre-Planckian) epoch took place in the beginning of cosmological expansion, 
consequent quantum theory of gravitation is necessary to analyze the PCS 1 . But such theory 
does not exist at present. By supposing that the Planckian epoch was absent by evolution of 
the Universe, classical theory of gravitation (with quantum description of gravitating mat- 
1 Recently the PCS was discussed in the frame of string theory (see [3, 4] and refs given here). 
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ter)can be used by analysis of the PCS. There were many such attempts to resolve the PCS 
corresponding to the following two approches (see [5] and refs given here). 

In the frame of the first approach one supposes that energy dominance conditions are valid 
for gravitating matter with extremely high energy densities and pressures at the beginning 
of cosmological expansion. In this case GR must be replaced by certain other gravitational 
theory by description of gravitating matter at such extreme states. Gauge theories of gravi- 
tation (GTG) and at first of all the Poincare GTG are an important generalization of GR. 
Cosmological equations for homogeneous isotropic models deduced in the frame of GTG 
lead to regular in metrics cosmological solutions in the case of usual gravitating matter [6]. 
These solutions differ from corresponding solutions of GR only at extreme states of gravi- 
tating matter near certain limiting energy density, where gravitational interaction has the 
character of repulsion instead of attraction. 

The second approach to resolve PCS concerns the analysis of description of gravitating 
matter at extreme conditions and examination of energy dominance conditions for corre- 
sponding states 2 . Unified gauge theories of strong and electroweak interactions with sponta- 
neous symmetry breaking are the modern base of matter description at extreme conditions. 
Inflationary cosmology as important part of the theory of early Universe was built by using 
gauge theories of elementary particles [8,9]. A number of problems of standard Friedmann 
cosmology were resolved in the frame of inflationary cosmology. Most inflationary cosmo- 
logical models discussed in literature are singular and their study is given from Planckian 
time. The radical idea of quantum birth of the Universe was introduced in order to avoid the 
PSC. Note that during quasi-de-Sitter inflationary stage gravitating matter in inflationary 
models is approximately in the state of so-called gravitating vacuum, for which pressure p 
and energy density p > are connected in the following way p = — p and energy dominance 
conditions are not valid. As it was shown in Refs. [10, 11] gravitating vacuum with suffi- 
ciently large energy density can lead to the vacuum gravitational repulsion effect (VGRE) 

in the case of systems including also usual gravitating matter, that allows to build regular 
2 Such analysis is important also in connection with the problem of dark matter and dark energy [7]. 
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inflationary cosmological models. At first the VGRE was discussed in the frame of Poincare 
GTG [10] in the case of homogeneous isotropic models including radiation and gravitating 
vacuum with p = const. Because cosmological equations of Poincare GTG are valid in the 
frame of the most general GTG — metric-affine GTG, the VGRE can take place also in 
metric-affine GTG [12, 13]. Regularizing role of gravitating vacuum in GTG was analyzed 
in Refs [14-17]. The VGRE in GTG can lead to a bouncing solutions in the case of closed, 
flat and open models. It is important that the VGRE can lead to the bounce also in GR 
in the case of closed models [11, 5] (see also [26] and refs given here). Regular bouncing 
closed models filled by linear massive scalar field were studied in the frame of GR in Refs. 
[18, 19]. By using values of some parameters of the Universe, the probability of regular 
solutions for such models was estimated to be very small [18] 3 . This estimation is not valid 
for inflationary models discussed below, for which scalar fields dominate during very small 
time intervals. 

By using some simplest potentials for scalar fields regular inflationary cosmological mod- 
els were discussed in the frame of GR as well as GTG in Refs. [22-25] . It was shown that 
limiting energy density and temperature at the bounce in inflationary models can be es- 
sentially smaller than the Planckian ones. Obtained regular inflationary solutions contain 
the stage of transition from compression to expansion and quasi-de-Sitter inflationary stage. 
After inflation scalar fields oscillate and are transformed into elementary particles, as result 
the transition to Friedmann regime takes place. Formal integration of cosmological equa- 
tions leads also to quasi-de-Sitter stage of compression [22-25]. However, this stage in GR 
is unstable and small change of variables at compression stage can lead to singular solution. 
This means that the building of regular inflationary models is possible, if physical conditions 
leading to the bounce take place at the end of cosmological compression (see below). 

By using results of numerical investigation performed in Ref.[27], the analysis of regular 
3 At the first time cosmological solution with initial de Sitter state limited in the past was discussed in Ref. 
[20]. The problems - what origin has initial vacuum state, what was before this state - were not considered 
in [20]. The cosmological solution with initial vacuum state not limited in the time in the past, proposed in 
Rcf.[21] does not lead to resolution of PCS in the future. 
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inflationary models in GR and GTG in dependence on conditions at a bounce is given below 
in present paper. 



2 Regular inflationary cosmological models in GR 

Let us consider homogeneous isotropic models filled by scalar field minimally coupled with 
gravitation and usual matter. We suppose that interaction between them is negligibly small. 
Then energy density p and pressure p can be written in the following form 

P = \<p 2 + V(<j ) )+p m , P=\<p 2 -V(<f ) )+p m , (1) 

where a dot denotes differentiation with respect to time, index m corresponds to usual 
matter, and V(<f>) is scalar field potential. Evolution of considered models in GR is described 
by Friedmann cosmological equations 

" + " 2 = 3^5 ( W) ~ ^ ~ \ (Pm + 3Pm) ) ' (3) 
where k = +1,0,-1 for closed, flat and open models respectively, R is the scale factor, 
H = R/R is Hubble parameter, M p is Planckian mass. (The system of units with h—c—1 
is used.) By virtue of negligible interaction of scalar field with usual matter the conservation 
law 

p + 3H(p + p) = (4) 

leads to equation for scalar field 

* + !Utf~V (V = Q (5) 

and 

p m + 3H(p m + p m ) = 0. (6) 
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In the case of closed models (k=+l) Eqs. (2)-(3) can lead to regular bouncing solutions. 
Putting t — at bounce (H(0) = 0, H = H(0) > 0), we see bouncing solutions take place, 
if the value of scalar field potential Vo = V ((f)(0)) (which plays the role of gravitating vacuum 
energy density) is sufficiently large and satisfies the following inequality 

v o > ^0 + \ (PmO + 3p m0 ) (7) 

(index corresponds to t = 0). For definiteness we will consider later the matter in the 
form of ultrarelativistic matter or radiation (p m = p r , p m = |p r , index r corresponds to 
radiation), then Eq. (3) takes the following form 

and according to (6) we have 

p r R A = const. (8) 

To obtain a bouncing solution we have to take some initial conditions for scalar field (4>o, 
0o) an d radiation p r0 satisfying according to (7) the following relation 

V - 4>l - PrO > (7') 

Admissible values of (fio, <p and p r0 must satisfy restrictions: V < \4>o ~ ^-^Ipi Pro ^ 

IMp , by which quantum gravitational effects are not essential and our classical consideration 
is valid. Minimum value of R of scale factor at bounce we obtain from Eq. (2) 



Ro 



71 (it 1 )2 

V + o^0 + PrO 



3M p 2 V 2 



(9) 



P o-^ 4 

By using initial conditions at bounce we integrate Eqs. (3') and (5) putting p r = ' in 

R 

accordance with (8). 

Now let us discuss the most important features of regular inflationary models in GR, 
by taking into account numerical results obtained in the case of the simplest scalar field 
potentials V 1 = ±A0 4 (A = 10" 14 ) and V 2 = \m 2 <\> 2 in Ref. [25]. 
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a) Transition stage from compression to expansion 



At first we will consider bouncing solutions in the case when the value of rj = Vq — 4% — 
p rQ > is not near to zero. Then as numerical analysis shows, the Hubble parameter H(t) 
and scalar field (pit) (0 O 7^ 0) vary in linear way during transition stage from compression 
to expansion, and derivatives H and <f> decrease rapidly and tend to zero at the end of this 
stage (t ~ ti), moreover radiation energy density p r — > at t ~ t\ (see Fig. I). 4 According 

10" 3 R 10 3 H 




Figure 1: Dynamics of functions R(t), H(t), <f)(t) and pit) during transition stage from 
compression to expansion in the case of initial conditions: V = lCr 6 Mp (0 O = 141.42M p ), 

00 = y/5W Q , p r o = O.Wq 



to Eq. (3') we have 

H(ti) = J^-, Vfa) (0! = ^0) (10) 

4 A11 figures in this paper are given in system of units h = c = M p = 1 for the potential V — |A</> 4 
(A = 10- 14 ). 
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and the change of Hubble parameter during time interval Ati — ti — t (t — 0, H — 0) is 



AH X = ^ = H Ah 



8tt 
3M$ 



V 



00 - Pro) 



Afi. 



(11) 



Because 0i = O + <f>oAti and hence according to (10) Hi = ,1 ^ V I 1 + 1 ^° Ati ) we 

3M~ \ 2 Vq 



find from (11) 



V 



(12) 



V ~ 00 _ PrO 



2v% 

Analogously we can determine time interval At 2 of transition stage before a bounce 

% 



AU 



57T 



3M2 



\ IV- 

V- 00 -PrO J +-^=00 



2x/V 



Then duration of transition stage is 



At tr = Ati + At 2 = 



(Mf» (^-0^ "Pro) 



8tt 
3Ml 



Vr 



— 00 — Pro) 



2 IV' 2 . 

--—0o 
4 V ° 



(13) 



and value of the Hubble parameter at the end of transition stage (at the beginning of 
inflationary stage) is 



8tt 



Vn 



- 00 - Pro) 



^ (~ 12 \ 1 V 1 



(14) 



3M p 2 



V - 00 - PrO - 



< 4 / (V - 0o - Pro) , we have 



3M 2 



(15) 



The value of Hi (15) depends on O only and does not depend on O and p r o- 
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By decreasing of value of rj, the linear time dependence of H(t) is violated and duration 
of transition stage increases (see Fig. 2). Formulas (13)-(14) are not applicable in this case, 
but expression (10) for the Hubble parameter at the end of transition stage is valid. 



10 J H 




Figure 2: Dynamics of function H{t) during transi- 
tion stage from compression to expansion in the case 
of initial conditions: Vo = 10~ 6 M 4 , O = VO-5Vo, 
p r o = 0A9V 

In the case of very small values of rj — > solutions change their character — inflationary 
stage and oscillations of scalar field vanish, solutions become unstable. For example, in the 
case of potential V = jXcj) 4 (A = 10~ 14 ) we have inflationary solutions for initial conditions: 
a) V = 1(T 14 M 4 , O = 0, p r0 = 0; b) V = 0.1M 4 , O = 0, p r0 = (l-l(r 6 )\/ ; c) V = 0.1M 4 , 
0o = a/(1 — 10 _3 )Vo, Pro — 0. However, we don't have inflationary solutions by the following 
change of these conditions: a) V = 1(T 15 M 4 (0 O = 0.8M P ), O = 0, p r0 = 0; b) V = 0.1M 4 , 
O = 0, p^ = (1 - 10- 7 )V ; c) Vo = 0.1M 4 , O = ^(1 - 10^)Vo, p r0 = 0. 

b) Inflationary stage 

During quasi-de-Sitter inflationary stage we have H <^ H 2 , <j) 2 <^ V ((/>), <f> <S V , p r ~ 
and term R~ 2 in Eq. (2) is negligibly small. Then according to Eq. (3') (and Eq. (2)) the 
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Hubble parameter at inflationary stage is equal to 



»- d -T- ik vw (16) 



and Eq. (5) for scalar field takes form 



Values of scalar field 0! = <j)(ti) and scale factor R x = R(ti) play the role of initial condi- 
tions by integration of Eqs. (16)— (17) for inflationary stage. For given potential V(4>) the 
integration of Eq. (17) leads to the function <j)(t). In the case of potential V\ = \\(p A we 
obtain 

0(t) = ie V 6 - pK 1J (18) 
and in the case of potential V 2 = \m 2 4> 2 we have 



= ^= "»(*-*!)■ (19) 



By using explicit form of scalar field 4> we can find the scale factor by integrating Eq. (16) 
In the case of functions (18)— (19) we obtain 

4tt 



R(t) = R x exp 



nMl 



(20) 



where n = 4 and n = 2 in the case of potentials V\ and Vi respectively. Expressions (18)-(20) 
have the same form as in chaotic inflation [8] , but in considered theory the values 0i and R\ 
are determined by initial conditions at bounce and are not arbitrary. 

At the end of inflationary stage Eqs. (16)— (17) are not valid, and behaviour of scalar 
field and Hubble parameter can be determined from Eqs. (3') and (5). After inflation the 
Hubble parameter is small and decreases, and scalar field changes as damped oscillations 
(Fig. 3). Oscillations characteristics (amplitude, frequency) depend on potential parameters 
and practically do not depend on initial conditions at bounce. Such dependence appears in 
the case of bouncing solutions with small values of r\ near their stability boundary [27]. But 
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Figure 3: Dynamics of functions H{t) and <f>(t) after inflation (Vo = 10 6 M^, <p = y/0.3Vo, 

PrO = 0.1^). 

noted solutions do not describe inflationary models and are not interesting for inflationary 
cosmology. Note, that relative part of such solutions in comparison with bouncing infla- 
tionary solutions is small. As result approximately relative part of bouncing inflationary 
solutions on the plane (0 O , 0o) satisfying inequality (7') (if p r0 = 0, < V((f) ) < 1M 4 , 
< \4>l < IM*) is equal to: a)^ for potential V x = ±A0 4 , b) ^ for potential V 2 = \m 2 <p 2 
[25]. 



3 Regular inflationary cosmological models in GTG 



Cosmological equations for homogeneous isotropic models obtained in the frame of gauge 
theories of gravitation have the following form 



k 

B 2 



d_ 

dt 



In 



Ry/\l-(3(p-3p)\ 



7T P - f (p - 3p) z 



3M p 2 l-p( p -3p) ' 



R + R(hiy/\l-P(p-3p)\)'] 47r p + 3p+ !( p _ 3 ^ 



(21) 



(22) 



R 3M p 2 l-P(p-Zp) 

Here (3 is indefinite parameter with inverse dimension of energy density. At first Eqs. (21)- 
(22) were deduced in Poincare GTG [6], and later it was shown that Eqs. (21)-(22) take 
place also in metric-affine GTG [12, 13]. The conservation law has usual form (4). Be- 
sides cosmological equations (21)-(22) gravitational equations of GTG lead to the following 
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relation for torsion function S and nonmetricity function Q 

S-\Q = -\±\n\l-0(p-3p)\. (23) 

In Poincare GTG Q = and Eq. (23) determines the torsion function. In metric-amne 
GTG there are three kinds of models [13]: in the Riemann-Cartan space-time (Q = 0), in 
the Weyl space-time (S = 0), in the Weyl-Cartan space-time (S ^ 0, Q ^ 0, the function 
S is proportional to the function Q). Eqs. (21)-(22) coincide practically with Friedmann 
cosmological equations of GR and noneinsteinian characteristics S and Q are negligible if 
energy density is small \P(p — 3p)| <C I. The difference between GR and GTG can be 
essential at extremely high energy densities \/3(p — 3p)\ > l. 5 In the case of gravitating 
vacuum with constant energy density p v = const > cosmological Eqs. (21)-(22) are 
reduced to Friedmann cosmological equations of GR and S — Q — 0, this means that de 
Sitter solutions for metrics with vanishing torsion and nonmetricity are exact solutions of 
GTG [28,29] and hence inflationary models can be built in the frame of GTG [22-24]. In the 
case of values of (3 to be not large in system of units with h = c= M p = 1 ( f3(AV — 4> 2 ) <1) 
inflationary models of GTG practically coincide with that of GR [24]. It is interesting to 



P(4V - 2 ) 



»1). 



investigate inflationary models in GTG in the case of large values of (3 

In order to analyze inflationary cosmological models in GTG let us consider systems 

including scalar field and radiation, for which quantities p and p are given by (1), where 

Pm = Pr, Pm, = \pr- Then cosmological equations (21)-(22) by taking into account Eq. (5) 
5 Ultrarelativistic matter with equation of state p — is exceptional system because Eqs. (21)-(22) are 
identical to Friedmann cosmological equations of GR in this case. 



11 



can be written in the following form 



k 



Z 2 + {# [l - 2/3(2V + 2 ) - 3/3V'<j>} 



H 



1 - 2{3(2V + 2 
+ 12(3H4>V 



3M* 
Z + H 2 

1 - 2(3(2V + 2 
8tt 



Pr + \& + V-\p[ \Y 



| 1 - 4p{V - 40 2 ) Z-18/3 2 4 } 



Z, 



(24) 



3/5 (F 



/ 2 



2 T W2 



3M2 



1 



\/-0 2 -p r -^(4V-0 2 ) 2 



(25) 



where Z 



P(W 



. Relation (23) takes the form 



13^ ^ + V 



(26) 



1 - /3 (w - 2 ) ' 

Eqs.(24)-(25) lead to essential restrictions on scalar field by evolution. So in the case k = 
0,+l from Eq.(24) follows that Z > 0, if j3 < 0. Because at the bounce (t = 0) we have 
H(0) = 0, from (24) follows 
k 



in 



1-2 1 , 

Pro + 2^0 +y -^/3(4Vb 



(4^0 - 



(27) 



(Zo = 1 — /3(4Vo — 0o))- The relation (27) determines minimum value of scale factor Rq for 
closed and open models, and in the case of flat models gives the dependence between </> , </> 
and p r Q. From (25) the time derivative H = H(0) at the bounce is 

1 



+ 3/5 



(V"^l - V'l)Z + 6f3<j> 2 V'l\ } [l - 2/3(2Vo + 2 



(28) 



Then the condition H > determines permissible values of (fio, <fio, Pro at the bounce in 
dependence on indefinite parameter (3. In order to analyze the dependence H = Hq{($) we 
write expression (28) in the form 



H = a h h ( P 2 + -P + - ) (P - Pi)' ' ; ' M 
a a 



(29) 
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where (3\ = |(2Vo + 0g) 1 and P2 = (4Vo — 0g) 1 are particular points and 
6 = 3l^ ( 4K) " ^) ( 4Pr0 + 5 ^ " 8V °) " 3 " ' 



If 6 2 — Aaf > 0, the formula (29) takes the form 

Ho = afafo (p - /3< 0) ) (p - /3f ) {p - {(3 - ft)" 1 , 



where (3[°l = ^ — . In the case of large (in module) values of/? (\(3 (^4:V — 0oj | 3> 

1) we have H « a/?i/?2 and the bounce condition H > leads to the following relation 

2 " Uv -^) 2 + 3(v^-v yi) + ^^ x >o w) 



The relation (|29IJ) is compatible with the following condition 

4V - > 0, (30) 

which follows from discussed above restriction on the value of Z in the case large (in module) 
negative value of /3 (/3 < 0). Corresponding inflationary models are regular in metrics, Hubble 
parameter and torsion (nonmetricity) . (Note that according to Eq. (24) the right-hand part 
of Eq. (26) tends to \H at Z — > 0). Unlike the condition (7') of GR, inequality (30) does 
not include radiation energy density p r . In contrast to GR, the appearance of VGRE in 
GTG (with large negative (3) does not depend on radiation, although the contribution of 
ultrarelativistic matter to energy density of the Universe at the bounce can be essentially 
greater than contribution of scalar fields. 

Let us consider cosmological equations of GTG (24)-(26) in approximation of large neg- 



ative P 



(3 [AV - 2 ) > 1), by supposing that p r + \(j) 2 + V < \(3\ [aV - (j) 2 ^ 6 We have 



3 This assumption does not exclude that radiation energy density can dominate at the bounce p r 



y+2 



U2 



13 



R 2 



+ 



2H[2V + d) 2 ) + 3V 



2tt 
3M 2 



AV 



2 ), 



(31) 



H (W + 2 ) (AV - 2 ) + H 2 (V 2 - 34V<j) 2 - 4 ) - 12HV'<j> (w + 2 ) 

3 



+ \ (V'V 2 - V" 2 ) (4V - 2 ) - 9V" 2 2 = ^ (4V - 2 ) , (32) 



4 V 2 4V-0 2 



(33) 



Eqs. (31)-(33) do not include radiation energy density, which is not essential for dynamics 
of regular inflationary models in considered case. According to Eq. (31) the scale factor at 
the bounce Rq is determined from 

2 



2tt 
37Wf 



(4^o - <Pl) ■ 



(34) 



The value of Rq determined by (34) coincides with that of GR (formula (9)) only if p r0 = 
and 0o — 0. According to Eq. (32) we have 



Hn 







4V^-0 2 ) +9^ 2;:; 



'0 



(35) 



that corresponds to (29) in the case of large values of \(3\. In accordance with Eq. (34) for 
closed models the following inequality at a bounce takes place 

2?r / -\3 



3M 2 



(4VW 2 ,)' 



9V ' 2 <j) 2 > 0, 



and in the case of open and flat models we have 



2 " V o -0 2 ) 3 -9VV 2 2 <O. 



3M2 

Because for open and flat models the derivative does not vanish, their evolution is essen- 
tially asymmetric with respect to the bounce point. 
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At the beginning of the inflationary stage (t ~ t\) the Hubble parameter H(ti) from (32) 
is equal to 

mi) = ^mf^ + ^ m = nm) ■ m 

In the case of linear time dependence of Hubble parameter during stage of transition from 
compression to expansion we can simply estimate the time interval Ati by using relations 
(35)-(36) and then determine the values of scale factor Ri = Ro{l + H Ati) and scalar field 
0i — 0o + 4>oAti at the end of transition stage (as it was made for inflationary models in GR 
in Section 2). Numerical analysis of regular inflationary models in GTG with large \/3\ shows 
that the most important difference of such models in comparison with that of GR concerns 
the transition stage from compression to expansion and the stage after inflation. As was 
noted above the region of initial conditions at the bounce of regular inflationary solutions 
is essentially more wide in GTG than in GR (compare (7') with (JHJ), (30)). In Fig. 4 the 
evolution of different characteristics of regular inflationary model in GTG {(5 = — 10 20 M~ 4 ) 
is given by choosing of initial conditions, by which a bouncing solution in GR is impossible. 
From Fig. 4 we see that the Hubble parameter and torsion as well as scalar field oscillate 
after inflation. Amplitude and frequency of scalar field oscillations are smaller in comparison 
with that of GR (for given potential) and decrease by increasing of 

Conclusion 

Considered regular inflationary models can be interesting for inflationary cosmology, if the 
realization of physical conditions leading to a bounce can take place at the end of cosmological 
compression. In the frame of GR this means that the greatest part of energy density of the 
Universe has to be determined by scalar fields (see (7')). However, in the framework of GTG 
the realization of the bounce is possible, when usual ultrarelativistic matter dominates at 
the end of cosmological compression, and scalar fields satisfy more weak restriction (relation 
(30)) — its kinetic energy density can be greater than the scalar field potential. Unlike 
GR, in the frame of GTG the elimination of cosmological singularity and realization of the 
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Figure 4: Dynamics of characteristics of inflationary model in GTG with (3 = — 10 °M 1 
during transition stage and stage after inflation (V = 10~ 6 M^, O = y/2V , p r0 = lO 5 ^). 



bounce in inflationary models are ensured by cosmological equations (24)- (25). 



Acknowledgements 

I am grateful to my pupil Alexander Garkun and my son Andrey Minkevich for technical 
help by preparation of this paper. I am thankful to my wife Eleonora for her support. 



16 



References 

[1] R. Penrose, Structure of Space- Time. (Mir, Moscow, 1972). 

[2] S. Hawking, G. Ellis, The Large Structure of Space- Time. (Nauka, Moscow, 1979). 

[3] M.Gasperini, G.Veneziano, Phys. Rept 373, 1 (2003) ; \hep~th/020713(\ 

[4] H. Liu, G. Moore, N. SeibeTg, \gr~qc/0301001\ 

[5] A. V. Minkevich, Acta Phys. Polon. B29, 949 (1998). 

[6] A. V. Minkevich, Vestsi Akad. Nauk BSSR. Ser. fiz.-mat., no. 2, 87 (1980); Phys. Lett. 
A80, 232 (1980). 



[7] M. S. Turner, \a~s~tw- ph/0207297 



A. D. Linde, Physics of Elementary Particles and Inflationary Cosmology. (Nauka, 
Moscow, 1990). 

[9] E. W. Kolb and M. S. Turner, The Early Universe (Addison- Wesley, 1990). 
[10] A. V. Minkevich, Dokl. Akad. Nauk BSSR. 29, 998 (1985). 

[11] A. V. Minkevich, Vestsi Akad. Nauk Belarus. Ser. fiz.-mat., no. 5, 123 (1995); in Proc. 
of the Third Alexander Friedmann Intern. Seminar on Gravitation and Cosmology, St.- 
Petersburg, July 4-12, 1995, p. 12-22. 

[12] A. V. Minkevich, Dokl. Akad. Nauk Belarus. 37, 33 (1993). 

[13] A. V. Minkevich and A. S. Garkun, Class. Quantum Grav. 17, 3045 (2000). 

[14] A. V. Minkevich, N. H. Chuong, Vestsi Akad Nauk BSSR. Ser. fiz.-mat., no. 5, 63 (1989). 

[15] A. V. Minkevich, I. M. Nemenman, Class. Quantum Grav. 12, 1259 (1995). 

[16] A. V. Minkevich, A. S. Garkun, Gravitation & Cosmology. 5, 115 (1999). 

17 



[17] A. V. Minkevich, G. V. Vereshchagin, Gravitation & Cosmology. 8, N4(24), 305 (2002). 
[18] A. A. Starobinsky. Pis'ma v Astron. Zhurn. 4, 155 (1978). 

[19] V. A. Belinsky, L. P. Grishchuk, Ya. B. Zeldovich,I. M. Khalatnikov, Zhurn. Exp. Teor. 
Fiz. 89, 346 (1985). 

[20] E. B. Gliner, I. G. Dymnikova, Pis'ma Astron. J., 24, 235 (1975). 

[21] H. J. Blome and W. Priester, Astron. Astrophys. 250, 43 (1991). 

[22] A. V. Minkevich, in Nonlinear Phenomena in Complex Systems. Proceedings of the 
Tenth Annual Intern. Seminar NPCS-2001, Minsk, 2001, pp. 231-236. 

[23] A. V. Minkevich, G. V. Vereshchagin, ibid., p. 237-240. 

[24] A. V. Minkevich, Intern. J. Mod. Phys. A, 17, 4441 (2002). 

[25] A. V. Minkevich, A. A. Minkevich, in Nonlinear Phenomena in Complex Systems. Pro- 
ceedings of the Eleventh Annual Intern. Seminar NPCS-2002, Minsk, 2002, pp. 209-213. 



[26] G. Ellis, R. Maartens, gr-qc/0211082 



[27] A. V. Minkevich, A. A. Minkevich, to be published. 

[28] A. V. Minkevich, Phys. Lett. 95 A, 422 (1983). 

[29] A. V. Minkevich, Dokl. Akad. Nauk BSSR. 30, 311 (1986). 



18 



